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We present an in-depth analysis of the spectral and temporal behavior of a long almost uninter- 
rupted INTEGRAL observation of Vela X-1 in Nov/Dec 2003. In addition to an already high 
activity level, Vela X-1 exhibited several very intense flares with a maximum intensity of more 
than 5 Crab in the 20-40 keV band. Furthermore Vela X-1 exhibited several off states where 
the source became undetectable with ISGRI. We interpret flares and off states as being due to 
the strongly structured wind of the optical companion: when Vela X-1 encounters a cavity in the 
wind with strongly reduced density, the flux will drop, thus potentially triggering the onset of the 
propeller effect which inhibits further accretion, thus giving rise to the off states. The required 
drop in density to trigger the propeller effect in Vela X-1 is of the same order as predicted by 
theoretical papers for the densities in the OB star winds. The same structured wind can give rise 
to the giant flares when Vela X-1 encounters a dense blob in the wind. Further temporal analysis 
revealed that a short lived QPO with a period of 6800 sec is present. The part of the light curve 
during which the QPO is present is very close to the off states and just following a high intensity 
state, thus showing that all these phenomena are related. 
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1. Introduction 

Vela X-l is a high mass X-ray binary (HMXB) consisting of the super giant HD 77581 and a 
massive (1.9 M 0 [27]) neutron star in a 8.964 day orbit [38]. The optical companion has a mass 
of ~23M 0 and a radius of ~30 Rq [38]. The neutron star is deeply embedded in the dense stellar 
wind of the donor HD 77581 (M* = 4 x 10~ 6 M G yr _1 ) [24]. X-ray lines indicate that this wind is 
inhomogeneous with many dense clumps [26] embedded in a far thinner, highly ionized component 
[31]. 

The neutron star has a long spin period of ~283 s [28]. The evolution of the spin period 
is best described by a random walk as expected for a wind-accreting system [40]. Although the 
source exhibits strong pulse-to-pulse variations, a pulse-protile folded over several pulse periods 
shows remarkable stability [34], even over decades [29]. At energies below 5 keV, the pulse-protile 
consists of a complex five-peaked structure, which transforms at energies above 20keV into a 
simple double-peaked pulse-protile [34] where the two peaks are thought to be due to the two 
accreting magnetic poles of the neutron star. 

With an X-ray luminosity of ~4 x 10 36 erg s -1 , Vela X-l is a typical high mass X-ray binary. 
Previous observations have shown that the source is strongly variable with reductions to less than 
10% of its normal value [20, 19, 17, 14], while periods of increased activity have also been observed 
during which the flux increases within an hour to a multiple of the previous value, reaching peak 
flux levels close to 1 Crab [19, 10, 16]. In this respect, Vela X-l is similar to sources such as 
4U 1700—377 and 4U 1907+09, for which low luminosity states and flares have also been observed, 
as is rather typical for wind-accreting systems [7, 37, 13, 1 1]. Although Vela X-l is a well studied 
object, only observations by INTEGRAL revealed that the flares in Vela X-l can be brighter than 
previously anticipated [20, 35, 21]. 

* Speaker. 
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Figure 1: a) ISGRI light curve in the 20-40 keV energy band and b) hardness ratio. The labels indicate the 
revolution number (from [20]). 


2. Data 

INTEGRAL observed the Vela region continuously for five consecutive INTEGRAL revo- 
lutions from revolution 137 (JD 2452970.86) until revolution 141 (JD 2452970.86) resulting in 
approximately 1 Msec of data (see Fig. [I]). 

We used all available science windows (ScWs) to be able to derive a contiguous light curve 
with as few interruptions as possible. Since Vela X-l is a bright source, the OS A can detect the 
source and determine its flux level accurately even when the source is at high off-axis angles. 
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Time [MJD-52971] 

0.14 0.15 0.16 0.17 0.18 0.19 



Figure 2: Close up of the light curve of giant flare 1 with a time resolution of 20 s. The peak is reached at 
MJD 52971.18 with 923 counts s _1 (corresponding to 5.2 Crab) in the 20 to 40keV band (from [20]). 


Apart from Vela X-l, also 4U 0836—429, H 0918—5459, the Vela Pulsar, and two sources first 
reported by INTEGRAL [5, 32] are detected. Vela X-l is widely separated from the other sources 
such that contamination of the spectrum of Vela X- 1 is of no concern. Data from JEM-X have 
not been used in this analysis due to the far smaller field-of-view of JEM-X and since Vela X- 1 is 
off-center in the observed field (Vela X- 1 was only within the fully coded field-of-view of JEM-X 
for less than ten out of the ~550 individual pointings). SPI data were also not used as SPI can not 
provide light curves on the short time scales required here. 


3. Data analysis 

Vela X- 1 was found in a highly variable state during the observation. While Vela X- 1 is known 
to be a variable source [19, 9], the behavior found in this observation [20, 35] is indeed extreme. 

Most importantly, on 2003 November 28 (JD 245297 1.67), INTEGRAL observed an extremely 
bright flare (flare 1 ; see Fig. [T|). During the flare, the 20-40 keV count rate increased from a SCW 
averaged pre-flare value of ~55countss _1 (~300mCrab, or 1.6 x 10 9 erg cm 2 s~') by a factor 
of more than seven to 405 counts s _1 (2.3 Crab) within only 90 minutes. 

After the peak, the flare decayed quickly to an intensity level of < 1 Crab and within ~ 1 1 h to 
a SCW averaged count rate of ~35 counts s 1 (200mCrab), somewhat lower than before the onset 
of the flare (see Fig. [I]). 

In the following, three more flares (flares 2 to 4, see Table [l]) were observed. All three flares 
were shorter and less intense than flare 1 on a science window averaged basis, but still reached 
SCW averaged intensities close to 1 Crab. 
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Table 1: Overview of the observed flares. See Fig. 1 for the numbering of the flares. The time is the onset 
of the flare. To obtain the peak fluxes f pea |., a light curve with a time resolution of 20 s was used. 7, vse is the 
time from the onset of the flare to the peak, while 7^ ota i is the duration of the flare. 


Flare 

Time 

[MJD] 

Duration 

[s] 

^peak 

[Crab] 

Trise 

Ttotal 

Remarks 

1 

52971.15 

11 200 

5.2 

0.15 

giant flare, spectral softening 

2 

52975.34 

5 200 

2.6 

0.83 

no spectral change 

3 

52976.50 

1800 

5.3 

0.28 

giant flare, very short 

4 

52977.15 

12 900 

1.9 

0.13 

spectral softening 

5 

52980.31 

31400 

3.9 

0.63 

high intensity state, no spectral change 


On 2003 December 7 (JD 245981.10), another intense flare was observed (designated Flare 5, 
see Fig. [I]). Unlike flare 1, during which the brightness of the source increased rapidly, it took ~8 h 
for flare 5 to reach its SCW averaged maximum 20^10 keV flux of ~1.2Crab. The decay lasted 
~5h until Vela X-l reached its pre-flare count rate of ~35 counts s -1 (200mCrab in 20-40 keV). 
Although quite bright, flare 5 is therefore significantly less intense than giant flare 1 , and also far 
longer, i.e. it is a high intensity state. 

The analysis of a light curve with a 20 s time resolution showed that the source reached a 
similar peak count rate in flares 1 and3: ~920countss _1 (5.2Crab) in flare 1 and ~930 counts s" 1 
(5.3 Crab) in flare 3. Flare 3 on December 3 was therefore also a giant flare. Flowever, flare 3 was 
significantly shorter: the entire flare lasted less than 2000 s (see Table |T]). 

Extending the analysis to the non-flaring parts of the light curve, we detected a quasi-periodic 
oscillation (QPO), similar to other accreting X-ray pulsars. The short-lived QPO with a period of 
~6820 sec appears to be quite regular and inconsistent with pure stochastic behavior (see Fig. |3). 
Subsequent period searches on the corresponding data subset clearly detect the period. We note 
that the quasi-periodic modulation shown in Fig. |] is far stronger than and inconsistent with the 
NOMEX effect, which can cause intensity variations from SCW to SCW, but not within a given 
SCW. 

Furthermore, we observed several off states, during which no significant residual flux was 
detectable by ISGRI (Fig. |4). The onset of these off states occured very suddenly. The luminosity 
of the source simply drops below the detection limit of ISGRI. At the end of the off states, Vela X-l 
switches instantly on again and immediately resumes its normal intensity level. Ah off states 
occured within 12 h from MJD 52981.0 and MJD 52981.5. 


4. Spectral evolution during the flares 


Although the source was extremely bright during the flares, meaningful spectral fits could 
not be obtained, as the exposure time was too short. We therefore analyze the hardness ratio (see 
Fig. [l]b), which is defined as 


HR = 


H-S 
H + S 


(4.1) 


where H is the count rate in the hard band (40-60 keV) and S the count rate in the soft band 
(20-30 keV). While the hardness ratio remained constant throughout most of the observation at 
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Time [MJD-52980] 

1.0 1.1 1.2 1.3 1.4 



2003 December 8 

Figure 3: Closeup of the light curve (with a time resolution of 283 s to remove the pulsations) where the 
temporary QPO is present. Note that during the trough between 2 h and 3 h, and especially following the 
quasi-periodic modulation, the count rate decreased several times to zero for a short time (from [20]). 


Time [MJD-52981] 

0.34 0.35 0.36 0.37 0.38 



2003 December 8 


Figure 4: Close-up on off states 3 and 4 of Vela X-l during which the source becomes undetectable by 
ISGRI and then turns on again within one hour (from [20]). 
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Figure 5: Hardness Intensity Diagram of Vela X-l; data from the eclipses have been excluded. Data from 
flare 3 are not incluved as INTEGRAL was in an engineering mode such that no spectral information could 
be obtained. The data points from the flares are indicated by individual symbols, (from [20]). 

~ —0.735, the hardness ratio significantly changed with the onset of flare 1: it dropped to —0.82 
and during the flare to —0.85 (see Fig. |I]b). The same for flare 4: the hardness ratio dropped from 
—0.72 to —0.84, the same level as in giant flare 1, although flare 4 reached only a third of the peak 
flux of flare 1. During flares 2 and 5, however, the hardness ratio did not change. 

We then used a “Flardness intensity diagram” (HID) to study the spectral evolution of the 
source (see Fig. |;). Most of the data points are centered around the average values of intensity and 
hardness ratio. The only exception are the data points from the flares, which are above the general 
cluster of data points. Due to the spectral softening, the data points of flare 1 and flare 4 are shifted. 
The softening, however, did not evolve during the flares, but the flares are softer than the average 
spectrum from the beginning until the end. 

5. Discussion 
5.1 The flares 

Vela X- 1 has always been known to be highly variable with time and to show intensity varia- 
tions of up to a multiple or a fraction of the original intensity on all time scales. Although Vela X-l 
has exhibited extensive flaring activity in the past, giant flares (as flares 1 and 3) had not been seen 
before. 

The analysis of the hardness ratio shows that there seem to be two types of flares: the first type 
(flares 1 and 4) shows dramatic increases in the count rate, the onset of the flare is very sudden, 
and the spectrum softens during the flare. The second type is similar to a high intensity state: these 
flares are longer and the spectrum does not change. 


♦ 

♦ Flare 1 

■ Flare 2 

a Flare 4 

♦ 

• Flare 5 

♦ — 



-0.9 -0.8 -0.7 -0.6 

Hardness 


7 


Highly structured wind in Vela X-l 


Ingo Kreykenbohm 


The mechanism behind the flaring activity, however, is not fully understood. It has been shown 
that a temproary disk may form in Vela X-l [36]. The disk collapses and the material is accreted 
onto the neutron star giving rise to a short flare. These predicted Hares would last from 15 to 60 
minutes, similar to the short flaring activity of Vela X-l. Furthermore, wind accretion is a highly 
instable process by itself: the accretion wake trailing the neutron star contains filaments which also 
produce flares when being accreted [2]. The shock trailing the neutron star oscillates creating the 
“flip-flop instability” which then produces inflows that repeatedly change their direction [23, 1]. 
The timescales of 45 min match some of the observed behavior very well, but fails to explain the 
long flares. 

As the local density in a shocked wind varies by a factor of 100 [15], which can explain the 
flaring X-ray luminosity [26]. Dense clumps trapped in an otherwise thin and more homogeneous 
wind might be responsible for long flares [22], when the clumps are being accreted. Such a clump 
can feed the neutron star with a significantly higher M than usual over several hours. In summary, 
we conclude that the observed long flares are due to a strongly structured OB star wind, while when 
Vela X- 1 is less active, the OB star wind is less structured. 

5.2 The off states 

In a similar way to the flaring activity, the off states [14, 19, 20] where the source is below 
the detection limit are remarkable. After the off state observed by RXTE in 1996 [19], the source 
resumed its normal, pulsating behavior without any transition phase. The offstates reported here 
also occured without a transition phase [20]. 

The reasons for these off states and the sudden reappearance of pulsations are not understood. 
So far, several ideas to explain these phenomena have been proposed, however, none can fully 
explain the observed off states, since all of these ideas require a significantly longer timescale, e.g. 
clumps in the stellar wind [6] would have to have an unrealistically high optical depth to completely 
block hard X-rays and can not pass the line of sight in a few seconds to explain the sudden turn 
on/off behavior of the source. Therefore other mechanisms must be considered. 

The wind of OB super giants is inhomogeneous and clumpy [39, 2]. The density of the stellar 
wind can vary by several orders of magnitude [30]. Not only clumps, but also holes, i.e. regions 
of strongly reduced density are present in the wind: there, the density is lower than the average 
density of the wind by a factor of 10 3 [30]. If the neutron star enters these holes, M would then also 
decrease by a factor of ~ 10 3 and the X-ray luminosity would be reduced accordingly. Furthermore, 
the density fluctuations predicted by these models occur suddenly [see Fig. 1 in 30] similar to the 
onset of the off states (see Fig. [I]). If M drops due to these density variations in the wind, the Alfven 
radius will increase due to the reduced ram pressure of the infalling gas. Once it is larger than 
the co-rotation radius, accretion onto the neutron star is inhibited, i.e. the X-ray source basically 
switches off. This scenario is commonly known as the propeller effect [12]. Since the propeller 
effect depends on the amount of infalling material, the Alfven radius is not constant. This effect 
was observed in GX 1+4 [4]: in very low luminosity states, no pulsations were observable, while 
the source was strongly pulsating in high luminosity states. 

Since the strength of the magnetic field of Vela X-l is known from the observation of the 
cyclotron lines [18], the critical flux limit for Vela X-l for the onset of the propeller effect can 
be obtainted (after [4]): Fx,p r0 peiier,VeiaX-i ~ 1.1 X 10 12 erg cm 2 s . Compared with the typical 


Highly structured wind in Vela X-l 


Ingo Kreykenbohm 


bolometric flux of several times 1 0 9 erg cm“ 2 s _1 , this critical flux is lower by about three orders 
of magnitude. This flux limit matches very closely the predicted density variations in the stellar 
wind of IO 3 - 5 [39]. 

We therefore conclude that off states could be caused by a sudden drop in M that allows 
Vela X-l to enter the propeller regime. Intensity dips, however, arc longer, show a smooth transi- 
tion, and exhibit photoelectric absorption of more than 10 24 cm~ 2 . These dips arc readily explained 
by a dense blob in the wind passing through the line of sight. 

5.3 Connection with SFXTs 

The similarity between the flares and off states in Vela X-l and the behavior of Supergiant 
Fast X-ray Transients (SFXTs) [33] is striking. SFXTs arc high mass X-ray binaries that show very 
brief outbursts on timescales of hours or even only tens of minutes, and then remain undetectable at 
higher energies for months between outbursts [25]. SFXTs should be rather bright persistent objects 
[8] since the neutron star is deeply embedded in the dense stellar wind of the optical companion, 
however, the accretion is inhibited by the propeller effect. It has therefore been proposed that 
SFXTs harbor a magnetar [3] and that the extremely strong magnetic field ( B > 10 15 G) effectively 
inhibts accretion unless the ram pressure of the infalling gas is high enough such that accretion 
becomes possible for a short time. The giant Hares and off states of Vela X-l are therefore similar 
to these outbursts: in both cases the accretion of a dense blob of material causes the outburst or 
Hare, while reduced material infall causes the source to switch off. In summary, Vela X-l and 
SFXTs arc rather si mi lar objects, however, SFXTs arc usually in the off state, while Vela X-l is 
usually in a normal accretion mode. 
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